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We study the ground state properties of an atom with nuclear charge 



P^ ! Abstract 

o 

f^ ■ Z and N bosonic "electrons" in the presence of a homogeneous magnetic 

^D , field B. We investigate the mean field limit N ^ oo with N/Z fixed, and 

^^ I identify three different asymptotic regions, according to B <^ Z^, B r-^ Z^, 

i-Q ■ and B ^ Z^. In Region 1 standard Hartree theory is applicable. Region 

II . 3 is described by a one-dimensional functional, which is identical to the 

'^ I so-called Hyper-Strong functional introduced by Lieb, Solovej and Yngva- 

C^ ' son for atoms with fermionic electrons in the region B ^ Z^; i.e., for very 

strong magnetic fields the ground state properties of atoms are indepen- 
dent of statistics. For Region 2 we introduce a general magnetic Hartree 
functional, which is studied in detail. It is shown that in the special case 
of an atom it can be restricted to the subspace of zero angular momen- 
5^ I tum parallel to the magnetic field, which simplifies the theory considerably. 

The functional reproduces the energy and the one-particle reduced density 
matrix for the full A^-particle ground state to leading order in N, and it 
implies the description of the other regions as limiting cases. 
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1 Introduction 

The ground states of atoms with many electrons in magnetic fields have been 
studied in ||LSY94a|, |LSY94b|, PSY00||, and their energies have been evaluated, 



exactly to leading order, as some of the physical parameters tend to infinity. 
The atoms have been modeled by the nonrelativistic quantum mechanics of N 
fermionic electrons, with an unmovable pointlike nucleus of charge Z in a homo- 
geneous magnetic field of strength B. In order to shed some more light onto the 
interplay of the involved laws of physics, we investigate the effects of changing 
one of them: What would happen, if the electrons were bosons? 

So we study the ground state of the Hamiltonian - written in appropriate 
units - 
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where we set 

^B.i = HV,+5a(x,)f. (1.2) 

The vector potential is given by a(x) = e" x x/2, where x G M^, e" is the unit 
vector parallel to the magnetic field in ^-direction. This Hamiltonian acts on the 
symmetric subspace of £^(]R^^, d^^x). We subtract B for every particle because 
we are interested in the binding energy, which is now equal to the ground state 
energy E{N, Z, B) = inf spec Kn,z,b- 

In the study of asymptotics, as B and Z tend to infinity, we find a division 
into three different regions. They are - always in appropriate units - characterized 
by: B '^ Z^, B r^ Z^, B ^ Z^. This is in contrast to atoms with fermionic 
electrons, where five different regions have been identified: B <^ Z^^^, B ~ Z^^^, 
Z^l^ <^ B <^ Z^, B r^ Z^, B:^ Z^. See ||LSY94a|l and references therein. 



A simple heuristic argument: The length scale, which is typical for the quan- 
tum effects of a single particle in the magnetic field, is ~ 5"^/^. Typical energies 
are the differences of the Landau levels, 2B. On the other hand, the length scale 
typical for a particle in the Coulomb potential only is ~ Z^^ and hence the typ- 
ical energy range is ~ Z^. In Region 2, where B ~ Z^, the magnetic and the 
Coulombic effects are therefore of the same order of magnitude. In Region 1, 
where B '^ Z^, the Coulomb effects dominate in all directions. Magnetic effects 
will not contribute in leading order. In Region 3, where B ^ Z^, the magnetic 
effects dominate the dynamics perpendicular to the magnetic field. The electron 
with low energy is confined to the lowest Landau band, the typical wave func- 
tions are squeezed to needles with diameter ~ 5~^/^. (See [|AHS81| , |FW94|| for a 
detailed rigorous treatment and for citations concerning the history of this prob- 
lem.) Turning from the one-body system to the A^-body problem, we remark that 
Bose statistics has no effect on the size of the ground states, if the pair interac- 
tions are ignored. Moreover, the repulsion of the particles is of the same order of 
magnitude as the attraction of the nucleus, ii N ^ Z. So the scaling properties 
of the lengths and energies per particle remain the same, and the distinction of 
the three regions for many bosons is the same as for a single electron. 

We remark, that we can moreover identify a Region 4 with strong magnetic 
field, where the nuclear charge Z is fixed. In this region the asymptotics is, in 
leading order, independent of the statistics, as has been noted in its evaluation 
in IIBSYOOj] . It is to be described by the model of one- dimensional atoms with 
delta-function interactions. 

Scaling: In the following we will use the parameters A = N/Z, (3 = B/Z"^ 
besides A^. By scaling x -^ x/Z the operator Z~^H ^ z b '^^ unitarily equivalent 
to 

with ground state energy E{N, A, (3) = Z^^E_{N, Z, B). We are interested in the 
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limit A^ ^ oo of N~^E{N,X, f3) with A fixed. In Region 1 this hmit is coupled 
with /? — i> 0, and in Region 3 with /? — > oo, while /3 is fixed in Region 2. The 
asymptotics of the atomic structure and of its energy for large N in the three 
regions is modeled by energy functionals in generalized Hartree theory. 

In Region 1 standard Hartree theory is applicable. The energy functional, a 
functional of the density p, is 

^""[P] = y"|Vp^/^(x)|2d3x-|±p(x)rf3x + D[p,p], (1.4) 

where 

Dlp,p] = U^^£.d^y. (1.5) 

^ J 1-^ y| 

Its ground state energy is 

E^(A)= inf £^[p]. (1.6) 

p. J p=^ 



It is known, ||BL83||, that 



hm ^EiN,X,0) = \E^{\). (1.7) 

We will extend this result to 

1.1. THEOREM (Energy asymptotics for Region 1). If N ^ oo and 

(3 = l3{N) —>■ with A fixed, then 

hm lE(iV,A,/?) = lE^(A). (1.8) 

N~*oo I\ A 



In Region 2 Hartree theory has to be generalized. The basic idea remains: 
the electrons occupy the ground states of an effective one-particle Hamiltonian 
with a mean field potential which has to be determined by self-consistency. Now 
in the presence of a magnetic field, the ground state of the effective Hamiltonian 
may a priori be degenerate, so that the electrons can be distributed over a larger 
set of states. To take this into account, one has to consider in general one-particle 
density matrices T in this Magnetic Hartree Theory. The energy functional is 

£f''[T]=TT[{H^-P)T]- /^pr(x)d3x + D[pr,pr], (1.9) 

where pr(x) is the density defined by F, with J Pr = Tr[r]. We define the Hartree 

energy E^^{X,l3) as 

E^'^{\,p)= inf ^^"[r]. (1.10) 

r, Tr[rl=A f^ ' ' 
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This general form of magnetic Hartree theory is necessary, if some extra ex- 
ternal potential is added, or if one considers, e.g., molecules. For, in the presence 
of magnetic fields, some well known facts of ordinary quantum mechanics are no 
longer true: the ground state may be degenerate, and, in the case of an axially 
symmetric system, it can happen that the energy is not minimized by states with 
zero angular momentum [p^077| , |AHS7^| . But for the atom without perturbing 
forces, it turns out that the theory can be reduced to the consideration of pure 
states, rank one density matrices, with a zero angular momentum component par- 
allel to the magnetic field. In this case (|1.9|) simplifies to the Magnetic Hartree 
density functional 

^rtP] = j {\^^P\' + ^r'p -(3p- ^p) rf^x + D[p, p], (1.11) 

where r is the radial coordinate perpendicular to the magnetic field. This func- 
tional is the restriction of ( p.. 91 ) to density matrices of the form \y/p){y/p\- We 
will show in the next section that both S^^ and S^^ have the same ground state 
energy and density, so one could alternatively define E^^ as the infimum of S^^. 
The energy asymptotics for Region 2 are stated in the following theorem: 

1.2. THEOREM (Energy asymptotics for Region 2). If N ^ oo with A 
and (3 fixed, then 

hm lE(iV,A,/?) = ^E^^(A,/5). (1.12) 

In Region 3 the ground state of the atom is squeezed into a needle with 
diameter - in the scaled coordinates - ~ /3~^/^. The Coulomb interaction of the 
confined particles acts along the needle effectively like a one dimensional delta 
function, with coupling constant ~ ln/3. It thus dictates the typical extension 
~ (ln/3)~^ of the ground state wave function in the direction of the field, and the 
typical energies as ~ (ln/5)^. This effective reduction to a one dimensional system 
has been discussed in ||LSY94ai |JY96| , |BSYOO|| ; see also ||BRW99|| for related 



studies. In the appropriate scaling, z ~ {\n(3)x% with x" the coordinate in 
the direction parallel to the magnetic field, the theory is a Hartree theory for a 
one dimensional model. It is identical to the theory for Region 5 of fermionic 
electrons, which has been studied in ||LSY94a| , including an exact solution of the 



ground state problem. Its energy functional has been called the Hyper Strong 
Functional; it is 

8''^\p] = j{^/'\z)fdz-p{Q) + \j p{zfdz, (1.13) 

with its ground state energy defined as 

^"^(A) = inf £HS[p]. (1.14) 

P, /P=A 

We will prove 
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1.3. THEOREM (Energy asymptotics for Region 3). If N ^ oo and 

j3 = (3{N) -^ oo with A fixed, then 

li-^^^ = T^^'W- (1-15) 

The difference between bosonic atoms in our Region 3 to fermionic atoms 
in the fermion-Region 5 is in the condition of apphcabihty of HS-theory. The 
Pauh principle demands one needle for each electron, electrostatics makes them 
lying side by side. The fermionic atom is thus a bundle of A^ needles, with total 
diameter - in unsealed coordinates and with unsealed parameters - ~ N^/'^B~^/'^. 
The condition for validity of HS-Theory is that the diameter of the atom, in the 
directions of x-*-, perpendicular to the field, is much smaller than the characteristic 
length of Coulombic quantum effects, ~ Z^^. This condition is therefore B ^ 
NZ"^ for fermionic atoms. Bosonic electrons may all occupy the same needle. For 
them, the particle number does therefore not appear in the condition, which is 
now B > Z^. 

In the investigation of the limits we exploit four principles: 

(i) Restriction to independent particles, 

(ii) Spatial concentration near the center, 

(iii) Concentration in the lowest Landau band, and 

(iv) High field limit of the Coulomb-interaction. 

The principle (i) is fundamental for the validity of the Hartree theories: The 
atom with many interacting particles will be compared to models with indepen- 
dent particles in an effective mean field. 

The spatial concentration (ii) had not to be stressed in systems without magnetic 
fields. In the studies of fermionic electrons, [[LSY94a] , PSY00|| , it has been proven 



as a consequence of the superharmonicity of the repulsive interactions. We will 

also use this superharmonicity, but in a different way: It implies the vanishing of 

the parallel component of the angular momentum in the state which minimizes 

the Hartree energy. Since one of the consequences is the absence of an "angular 

momentum barrier" , it can also be viewed as a spatial concentration of the bound 

electrons. 

The principles (iii) and (iv) are effective in Region 3, in the limit (3 — > cxd. 

In the investigations of bosonic "electrons" we could probably have mimicked 
the procedure of |P!jSY94a||, with some changes due to the Bose statistics. Our 



procedure relies in fact heavily on the same methods, but we combine them in a 
new way: The principles (ii), (iii) and (iv) mentioned above are studied for single 
particles in effective mean fields. In the study of many particle systems, we begin 
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with the reduction to independent particles. But this has to be done in a subtle 
way, anticipating the limits which have to follow. We do this by extending the 
method of pSYOg. 



The Hartree theories will be discussed in Sect. 0, including the restriction 



to zero angular momentum in Sect. 2.5. The confinement to the lowest Landau 



band, relevant for Region 3, is treated in Sect. |3]^ and the following subsection 
In Sect. § Hartree theory is proven to be the limit of many particle quantum 
mechanics, as it is formulated in the Theorems |1.1| , 1.2 and |1.3| . Subsection O 



treats the restriction to the independent particle model. Finally, some results on 
the states and on "Bose condensation" are presented in Subsections ^^ and ^?7[ 



In the investigations of the limiting procedures, we are interested in the phys- 
ical dimensions of estimates and bounds, not about numerics. We use "C" or 
^^const." for all the numerical constants. 



2 Hartree theory 

2.1 Definitions and basic properties 

Definitions: The Hartree functional without a magnetic field, S^[p] in (|1.4]) , is 
defined for non-negative densities p(x) e £^(R^,(i^x) with the restriction that 
every component of Vp^' ^ (x) is an element of C^ (M^ , d^x) . Analogously, the 
functional for Region 3, £^^^[p] in (|1.13|) , is defined for non-negative densities 
p{z) e C'^(R,dz) with the restriction d{p^^'^{z))/dz E C'^(R,dz). The functional 
for Region 2, £^^^[r] in ( |1.9|) , is defined for density-matrices, non-negative trace 
class operators F acting on £^(R'^, d^x), with the restriction of a finite magnetic- 
kinetic energy: 

Tr[HpT] < oo (2.1) 

The associated density pr(x) can be defined in £^(]R'^, d^x) as a norm convergent 
sum of integrable densities ^,^WkPk, by diagonalizing the density matrix as 

r = J2^k\i^k)M, (2.2) 

k 

with normalized ipk, and pfc(x) = |'?/^fc(x)p. 

The conditions of finiteness of the kinetic energies imply the finiteness of 
potential energies, in all three regions: the attraction is bounded by the kinetic 
energy, because of the boundedness of the Coulomb potential (delta function 
potential in Region 3) relative to the operator of kinetic energy, which is proven 
with the stability of the hydrogen atom. Moreover, the repulsion is bounded by 
attraction, because for cr(x) and p(x) both non-negative elements of i2^(M'^, (i^x), 

2D[a,p]<XsupAy[p], (2.3) 
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with A = Jcr(x)(i'^x and Ay[p] = f , ^ | p(x)(i^x. The analogous inequahty for 
the HS theory is 



a (z) p{z)dz < X sup p{y), (2.4) 

y 

for a{z) and p(2;) both non-negative elements of C^{W,dz) and A = J a'{z)dz. 

In the variational principles which define the Hartree energies, the restric- 
tions IIpIIi = A and Tr[r] = A can be weakened to ||p||i < A and Tr[r] < A, 
because one can always "move some charge to infinity". This will be used for 
the proof of the existence of a minimizer for the energy, and this in return means 
that "inf " can be replaced by "min" . 

In the following discussion we will explicitly study the magnetic Hartree the- 
ory. All the results which do not refer to the dependence on (3 are also valid - 
in their essence of physical meaning, with some changes in the mathematics - 
for the standard and the HS theory. The proofs can be transfered, keeping their 
structure, but changing the mathematical expressions. This is an indication, that 
the physics behind the arguments is often the same. We point out, in particular, 
that the delta potential has the same scaling properties as the Coulomb potential. 

Introducing more non-negative parameters, we will study the extended energy 
functional 



(2.5) 



^A7c,"[r] = XTriiHf, - (3)r] - A / ^pr(x)ci=^x + aX'D[pr, pr], 

and its ground state energy 

E^f{X,(3,C,a) = ^ inf ^^£,7,,„[r]. (2.6) 

By scaling one verifies that the energies are related through 

E^^{X,(3,C,a) = ^-E^^{^X,^^(3). (2.7) 

Monotonicity, convexity and concavity properties. 

(i) Since the functional ^AACat"*^] ^^ decreasing in (, increasing in a and jointly 
linear in {(, a), the Hartree energy E^ is decreasing in (, increasing in a 
and jointly concave in {(, a). 

(ii) Since the functional ^x]f^(^^a{^]/X is increasing in A and jointly linear in {(, A), 
the energy per unit charge of the electron cloud, E^^/X, is increasing 
in A and jointly concave in (C, A). 
Because of (pT^), these properties are actually equivalent to (i). 
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(iii) Since moreover the functional £^^^[r] is convex in F (it is even strictly 
convex in pr, since D[p, p] is strictly convex in p), and since the sets {Tr[r] < 
A} are ordered by inclusion, the Hartree energy E^^{X,j3) is decreasing 
and convex in A. For, if F„ with Tr[F„] = A and T„ with Tr[T„] = A are 
minimizing sequences, it follows, that E^^{{X + A)/2,/?) < inf„£^^^[(F„ + 
T„)/2]. 



The convexity property justifies the definition of a critical charge Ac, which 
a priori might be infinity, indicating the maximal charge which the nucleus can 
bind: E^^ is strictly decreasing for A < Ac, and constant for A > Ac- 

Moreover, the convexity of E^^ in A and the concavity of E^^/X give upper 
and lower bounds on d'^E^^/dX'^ , guaranteeing the existence of the chemical 
potential 

as a continuous function of A. 

Comparison with the hydrogen atom: 

The extended functional (^3|) with a = is obviously related to the energy of 
the hydrogen atom described with the Pauli Hamiltonian, where the magnetic 
moment of the electron serves for the subtraction of (3 in the ground state energy: 

E^f{X,(3,C,0) = XE'^^\f3,O- (2.9) 

We have the following bounds: 

2.1. PROPOSITION (Energies of hydrogen as bounds). For A > and 

a> 

E^f{X,(3,C,c^) > XE'^y'iPX) (2.10) 

E^f{X,(3,C,a) < Ai?^y<i(/3,C-Aa/2). (2.11) 

Proof. The first inequality follows from the strict positivity of D[p, p]. The upper 
bound can be given by choosing the projection onto the ground state of the 
Hamiltonian Hf^ ~ iC ~ ^<^/2)/|x|, as a test density matrix F. Then we apply 
the bound to the repulsion by attraction (p.3|), together with the observation, 
that the nucleus has to be at the point of the maximum of the potential ^y[p]. 
(Otherwise the energy could be lowered by shifting p). D 



We remark that the bound (|2.11| ) is of no use for Xa ~ 2( or larger. In this case 



one can use the monotone decrease of the energy E^^ in A, and bound E^ by 
m.m\{XE^^'^{P, ( — Aa/2)}. We will use these bounds in Subsection ^l2| . 

As an obvious consequence of these bounds, using also the continuity of the 
hydrogen energy in (^, we add 
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2.2. Remark (The limit A — >0). In the limit A — » the energy per unit 
charge, E^^{X, P, (, a)/\, converges to the energy of the hydrogen atom, E^^'^{P, (). 



The bounds of Proposition ^7T| will actually be used in a simplified form: 
2.3. LEMMA (Simple bounds). 

^^"(A,/5) > -(l/4 + /3)A, (2.12) 

E^«(A,/5) < -(l/4)A(l-A/2)2, for A < 2. (2.13) 

Proof. Applying the diamagnetic inequality ||S76|| to the lower bound in proposi- 



tion ( |2.1|) , we get 

Hp-(3--^>mi spec {Ho - (3 ~ -^) = -] - f3. (2.14) 

|x| |x| 4 

To get the upper bound, we apply Lieb's inequality (Theorem A.l. in ||AHS78[| ) 

inf spec (Hp-P + V) < inf spec (Hq + V) (2.15) 

to the upper bound in proposition ( |2.1| ). D 

2.2 Minimizers 

2.4. THEOREM (Existence of a minimizer). For each P > and A > 
there is a minimizer F^ for S^^^ under the condition Tr[r^] < A, i.e. 



nMH/ 



£;^vm(A,/3) = 4^^"[r"]. (2.16) 

Proof. We follow closely the proof of the analogous theorems 2.2 and 4.3 in 



[LSY94a|| . Let r„ be a minimizing sequence for S^^ with Tr[r„] < A. First 

Tj— v-f-/^ r y^ t~i r Iv I— I _l ii~i V'\ /''^ 1 1 irt /'I r\ /~\ t~i \r\ j— »■« r/-\ l--v/^/^(-\iTi~i/~\ >■ -IVl JTl. I ^r^ V-^ /-'t 1 1 irt /'I r\ /~\ t~i \r\ j— »■« r/-\ ri v^ /~\ 



note that TT[HpTn\ is bounded above, because £^^ [F„] is bounded above, and 
since the other contributions to the energy are bounded relative to if/?. Now 
we show that the magnetic-kinetic energy is bounded below by the 3-norm of p: 
Using the diamagnetic inequality (the prerequisite for the final result in ||S76|| ) 

(2.17) 



(2.18) 



and the decomposition of F as in ( p.2|) , we get 



Tr[(-^V + /3a(x))' T]>Y,Wn j M^n\ 



Moreover, using the Cauchy-Schwarz inequality. 



iVprr 



= ^2Wn\i'n\^\i^n\ 
n 

< A\S^Wn\^n\A\y^Wn\V\^n\A 



(2.19) 
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Because Vpr = 2pp Vpp this gives 



1/3 



Tr[(-^V + /5a(x))^r] > | Ivp^'f > 3 (|)'^' ^p^^ , (2.20) 

where we have used the Sobolev inequahty in the last step. 

We can then conclude that the corresponding sequence p„ = pr„ is bounded 
in L^ n C^ and pn is bounded in Ti,^. Therefore, for each p G (1, 3], there exists a 
subsequence, again denoted by p„, that converges to some poo weakly in C^ nC^, 
and pointwise almost everywhere. It follows from weak convergence that poo > 
and J Poo < ^- From Fatou's lemma we infer that 

liminfL'[p„,p„] > DfpocPoo]- (2.21) 

n— >oo 

Moreover, since |x|~^ G C^^'^ + C^^^ for every e > 0, and choosing p the dual of 
3 + e, we see that 

lim /"^Pn= [ tKpoo. (2.22) 

n^co J |x| y |X| 

Since the r„'s are trace class operators on £^, we can pass to a subsequence 
such that for some F^ 

lim Tr[F„A] = Tr[F^A] (2.23) 

n— >oo 

for every compact operator A. (Here we used the Banach-Alaoglu Theorem and 
the fact that the dual of the compact operators is the trace class operators). In 
particular, we have 

F„ - F« (2.24) 

in the weak operator sense. It is clear that F^ > 0. Now let (/)j G C^(]R'^) be an 
orthonormal basis for C^. Again by Fatou's Lemma 

Tr[F^] = V(0,|F^|0,) < liminf Tr[F„] < A. (2.25) 



In the same way one shows that 

Tr[(/J^-/?)FH] = Y.^{Hp-pf/'ct>,\T''\{Hp-pf"ct>,) 

j 
< liminfTr[(/7^-/3)F„]. (2.26) 

n— >oo 

It remains to show that ppn = poo- We already mentioned that for some constant 
C 

TT[HpT^] < C (2.27) 

for all n. It follows from ( |2.24| ) that 

(1 + H^y/'T^{1 + Hpf'^ - (1 + Hpf^V'^il + Hpf'^ (2.28) 
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weakly on the dense set of C^ functions. Since the operators are bounded by 
(p:27D , (|;28D holds weakly in C^. 



Now consider some / G C^ acting as a multiplication operator on C^. It is 
easy to see that / is relatively compact with respect to —A, i. e. /(I — A)~^ is 
compact. In fact, it is Hilbert-Schmidt, because the trace of its square is given 
by 

y"/(x)/(y)F(x-y)2rf3xrfV (2.29) 

with the Yukawa-Potential YijK.) = (47r|x|)~^exp(— |x|), and this is bounded by 



Young's inequality. From ||AHS78|| , Thm. 2.6, we infer that 

g = {l + Hpr'/'f{l + Hpr'/' (2.30) 

is compact (it is even Hilbert-Schmidt). Thus there exists a sequence gi of finite- 
rank operators which approximates g in norm. We have 



/ 



(p„ - Prn)/ 



= |Tr[(r„-r«)/]| 

< |Tr[(l + HpY/\T^ - rH)(l + H^y/'g.] \ 

+2(C+1)A||(7-(7,||, (2.31) 



where we have used (|2.27| ). Hence p„ -^ p^a in the sense of distributions. Be- 



cause we already know that p„ converges to p^o pointwise almost everywhere, 

we conclude that p-pn = poo- We have thus shown that there exists a F^ with 

Tr[F^] < A and ^^"[F^] < lim inf „^oo ^«^^ [r„] , from which we conclude that 
rMHrrHi _ zt'MHm ,Q^ 

^/3 



g^^[T'']=E^'^{X,p). n 



2.5. Remark. If F^ is unique, and given any minimizing sequence for Sg^^, the 
whole sequence converges weakly to F^. 

Although we cannot make an assertion about the uniqueness of F^ yet, we 
can state the 

2.6. PROPOSITION (Uniqueness of the density). The density p^ corre- 
sponding to the minimizer is unique. 

Proof. This follows immediately from the strict convexity of D[p,p]. D 

For A < Ac, the energy E^^ is a strictly decreasing function of A. So the 
minimizers for different A < Ac are different, and they have the normalization 
Tr[F] = A. No part of the electron cloud has to be moved to infinity. The strict 
convexity of D[p, p] implies now, for A < Ac, a strict convexity of E^^ as a 
function of A. 
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2.3 Some physical quantities and their interrelations 

Since the densities of the minimizers are unique, the contributions to the Hartree 
energy, Jp/|x| and D[p, p], are fixed. We denote them as attraction A and 
repulsion R, suppressing the dependence on the parameters. As a consequence, 
also K, the kinetic-magnetic energy, is fixed. Inserting a minimizer into (|1.9| ) we 
get 

E = K -A + R. (2.32) 

To deduce an analogue to the Feynman-Hellman theorem, we observe the fol- 
lowing inequality: Consider two different parameters, ( and (, with corresponding 
minimizers F and F, and their densities p and p. All the other parameters are 
fixed. Insert F into the functional ( |2.5| ) with parameter (, to conclude 



cMH_ 



[r] - (C - C) / ^PrH 



1 



^exf(A,/3,C,«)-(C-C) 



rP- 



(2.33) 



Together with the same argument, where ( and ( are exchanged, we get, for 

oc, 



1 E^fiX,f3,C,a) 



iP< 



E^^{X,(3,Ca) 



c-c 



< - 



fP- 



(2.34) 



We conclude, using the continuity (|2.22| ), where the p„ are now the unique den- 
sities of the minimizers for (n -^ C, that the Hartree energy is differentiable in (, 
and 



dE, 



MH 



dC 



-A. 



(2.35) 



C=i 



In the same way, using the inequality ( |2.21|) instead of ( p.22p , we show that 



dE. 



MH 
ext 



da 



R. 



(2.36) 



a=l 



Having established the differentiability in A, C, and in a, the scaling relation (|2.7| ) 
implies differentiability in (3. So, the magnetic moment 9, 



dE 



MH 



dj3 



(2.37) 



is well defined as a function of (3. 
Taking the partial derivatives of 

R 

-A 



in a and ^ at a = 1 and C = 1 gives 



-E + A/i, 
3E-\p- 2(36. 



(2.38) 
(2.39) 
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The relation ( p.38| ) is, because of the negativity of the chemical potential /i for 
A < Ac , a virial inequality, 

R<\E\, (2.40) 

turning into an equality at A = Ac- Inserting ( ^.38| ) into ( |2.39| ), we get 

K - \E\ 
(36 = -^. (2.41) 



The equation ( 2.38|) and the inequality ( 2.40 ) are also valid in the Regions 1 and 



3, where there is a scaling relation analogous to (|2.7|) . But this relation is without 
a (3, so instead of (p.41| ), in Regions 1 and 3 the well known virial equality 



K = \E\ 



A-R 



(2.42) 



holds (see also [|LSY944 |B84|| ). Moreover, for A = Ac we have \E\ : K : A : R 
1:1:3:1. 



2.4 The linearized theory 

Since the density p^ of a minimizer is unique and integrable, we can define a 
linearized Hartree functional by 



^,Tn[r] = Tr[i/«r], 

with the one-particle Hartree Hamiltonian 
Here the Hartree potential $^(x) is given by 

$H(^) = 1 _ H 1 _ 

X X 



(2.43) 



(2.44) 



(2.45) 



Note that H^ depends on A only via p^. 

Since $^ is in £^ + £^, it is relative compact with respect to Hp ||AHS78|| , so we 

know that the essential spectrum of H^ is [0, cxd). 



MH 



2.7. LEMMA (Linear Hartree functional). Let T^ be a minimizer of S^ 
under the constraint Tr[r] < A. Then F^ also minimizes S^^^ (under the same 
constraint) . 



Proof. (We proceed essentially as in ||LSY94a| ]). For any F 

Sf^[T] = £^l[T] + D[pr- p^ pr - P^] - D[p^, p«] 



(2.46) 
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Especially, for all S > 0, 

+<52D[pr-p«,Pr-p^]. (2.47) 

Now if there exists a Fq with Tr[ro] < A and ^^i^[ro] < ^^^ir^] we can choose 
6 small enough to conclude that 

£^[(1 - ^)r" + ^ro] < ^,^i^jr«] - D[p^, p«] = sf^n (2.48) 

which contradicts the fact that F^ minimizes S¥^. □ 

The Hartree equation, the Euler-Lagrange equation corresponding to the min- 
imization of the functional S^^^, is 

H^T = pF. (2.49) 

The ground state energy of H^ is given by the chemical potential fi{X, j3) = 
{E + R)/X. For the overcritical values of A > Ac the ground state energy p of H^ 
is 0, and there is just one density corresponding to the minimizer. 

We can now ensure that Ac = XdP) is not too small. In fact we state 

2.8. LEMMA (Lower bound on the critical charge). 

Ac(/?) > 1 for all (3>0. (2.50) 

Proof. For /3 = this was shown in | BL83 |, and Ac(0) was computed numerically 



to be 1.21 in B84 



Fix /3 > 0. We assume that Ap = Jp^ < 1, and will show that H^ has an 
eigenvalue strictly below zero. Using ip{r,z) = exp(— /3r^/4 — a\z\) with some 
a > as a (not normalized) test function we compute 

(^|ifH^) = ^a - 27r /" <l>"(r, z)e'^f''"-^^^'^rdrdz 
P J 

< ^a-y"p«(x)(0(O)-0(x))rf3x, (2.51) 

where 0(x) is the potential generated by the charge distribution exp(— |/5r^ — 
2a|2;|), i.e. 

r p-ki^-'^-MA 

0(y) = / —, ^d^^. (2.52) 

J |x-y| 

Note that 0(0) > 0(y) for y 7^ 0, so we can choose a small enough to conclude 
that H^ has an eigenvalue strictly below zero and binds more charge than A^. D 



Lemma 2.8 shows that in Hartree theory there are always negative ions. The 



following lemma gives an upper bound on the maximal negative ionization. 
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2.9. LEMMA (Upper bound on the critical charge). For some constant 
C independent of (3 

Ac(/5)<2 + lmin{l + AC(l + [ln/5]2)}. (2.53) 



Proof. This follows from Lemma ^^ and the bound on N^ given in [pOO|l (see 



Remark |4.10| below). D 



2.5 The minimizer has L" = 

Let H be the operator 

H = Hp-l3-^-K), (2.54) 

where the potential is given by 

$(x) = — - — *p. (2.55) 

The function p is assumed to be axially symmetric, non-negative, and p G £^n£'^, 
where q > 3/2. This implies that |x|~^*p is a bounded, continuous function going 
to zero at infinity. The operator H is essentially self-adjoint on C^(M^). Since 
$ G £^ + Cf, it is a relatively compact perturbation of if/3 ||AHS78|| , so we know 
that the essential spectrum of H is given by [0, cxd). The symmetry of p implies 
that H is axially symmetric, i.e. it commutes with the rotations generated by 
the parallel component of the angular momentum. 

For 771 > let now \E'm be a ground state of H, if there is one, with angular 
momentum L^^^^n = —fn'^m- Note that we can restrict ourselves to considering 
non-negative tti's, since (H \ L" = —m) is antiunitarily equivalent to {H + 2Pm \ 
L'I = m) by complex conjugation, so in the ground state m is certainly non- 
negative. Writing 

vl/„(x)=e-^"^^r"^/(r,z), (2.57) 

where {r,ip) denote polar coordinates for {x,y), we see that / is a ground state 
for 

Hm = -TT^ TT - TT^ + ^^ - (^ + 1)/^ - Hr, z) (2.58) 

or^ r or oz'^ 4 

on £^(]R^,r^'^(ix). If H^ has a ground state, it is unique and strictly positive 
( ||RS78|| ; to apply the theorems therein note that the first three terms in H^n are 



just the radial part of the Laplacian in 2m + 3 dimensions). So \E'm is the only 
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ground state of H with L" = —m. Moreover, / is a bounded, continuous function 



LL97|| , and hence 

|^m(r,0,2)| <const.r'", (2.59) 

for small r and some constant independent of z. In particular, \l/m(0) = if 
m > 0. 

We now are able to prove the main result of this subsection: 

2.10. THEOREM (m = in the ground state of H). Let H be a Harmlto- 
nian as in \2. 5^ ), with p axially symmetric and non-negative. If inf spec H is 



an eigenvalue, the corresponding eigenvector is unique and has zero angular mo- 
mentum, for all values of f3 >0. 



Proof. For p = 0, i.e. the hydrogen atom, this was shown in ||AHS81|| (and also 
in ||GS95|| ). So we can restrict ourselves to considering the case j p > 0. 



Let \l'm be a normalized ground state for H, with angular momentum —m. 
Define 

/(b) = -j |M/„(x)|2$(x - b)rf3x. (2.60) 

The function / is continuous and bounded. It achieves its minimal value at b = 
0, because otherwise one could lower the energy by translating |\l/mp. Strictly 
speaking, 

/(b) - /(O) = {^rn\H\^m) " (^,n|i^|^m) > 0, (2.61) 

with b = (pi, 1)2, h-^) and 

$„(x) = e-t^(''^"-''i^)vl/„(x + b). (2.62) 



The phase in ( p.62|) is chosen such that the kinetic energy remains invariant. 
(Note that translating Hp is equivalent to changing the gauge of the magnetic 
potential.) 

In the sense of distributions, 

A/(b) = 47r (|^^(b)|2 - |^„|2 * p_(b)) , (2.63) 

with P-(x) = p(— x). The function |^mP * P~ is pointwise strictly positive and 
continuous. Assume now that m > 0. Since, by ( |2.59| ), |\l/m(x)| < C|x|™ for 
some C > 0, there is an i? > such that 

A/(b) < for |b| < i?, (2.64) 

i.e. / is superharmonic in some open region containing 0. This contradicts the 
fact that / achieves its minimum at b = 0. As a consequence, a ground state of 
H must have angular momentum m = 0. D 

An analogous result holds also for the restriction of H to the lowest Landau 
band. This fact will be used in the proof of Theorem E^. 
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2.11. COROLLARY (m = in the ground state of UoHUo). Let Uq be 

the projector onto the lowest Landau band. If UqHUq has a ground state, the 
ground state wave function \I/o is given by 




^o(x) = y f^ exp [-^r'j ^Piz) (2.65) 

for some iplz) with J \ilj{z)\'^dz = H^E'olli- 

Proof. Mimicking the proof of the last theorem, we see that if UqHIIq has a 
ground state, the corresponding wave function \l/o = Hq^I/o has L" = 0. Hence it 
is given by ( p.65| ). D 

Let now p be the Hartree density p^, the unique density of the minimizer of 
£p^^ (depending on A and /3). Applying Lemma (|2.7|) and the theorem above, we 
immediately get 

2.12. COROLLARY (Uniqueness off"). The functional Sf^ has a unique 
minimizer F^ under the condition Tr[r] < A, which is proportional to the projec- 
tion onto the positive function \/p^. In particular, T^ has rank 1. Moreover, p^ 
minimizes the is density functional ( \1.1^ ) under the condition f p < \. It is C°° 
away from the origin, continuous ai x = 0, and strictly positive. 

The properties of p^ stated above follow from the fact that p^ minimizes S^^^, 
and therefore satisfies the variational equation 

-A^+((^r^-P-<^A^ = p^, (2.66) 



where p is the chemical potential of the MH theory. Note that Corollary |2.12 



proves the assertion made in the Introduction that both S^^ and S^^ have the 
same ground state energy and density. 

3 The limits of Region 2 

3.1 The limit of very weak magnetic fields 

The limit /3 ^ is rather easy to handle: 

3.1. THEOREM (Hartree energy for small p). In the limit P ^0, 

E^'^{\(3) = E'^{\)-(3\ + 0{(3^). (3.1) 
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Proof. Fix A > 0, and let Pj3 be the miniinizer of the density functional Sp^^ 
under the constraint j p < X. Using po as a trial function we get 



E^»(A, P) < Spipo] = i?^^(A, 0)-p\+^ I r'po. 



(3.2) 



It is well known that po falls off more quickly that the inverse of any polynomial, 
so / r^po is finite pBL81| . 



For the converse, we estimate 



E^«(A, /5) = Sf^'ipp] > S^iPp] - /3A > E^'^iX, 0) - /3X. (3.3) 



The observation E^"(A,0) = ^^(A) then leads to (O 



D 



3.2 Lowest Landau band confinement in Hartree theory 

We now show that for large f3 most of the charge is confined to the lowest Landau 
band. To do this, we write the Hartree functional as 



Sf^[T] = Tr 



H. 



P-n)^ 



|Tr2 



r®r- 



1 



(3.4) 



where Tr2 means the trace over the doubled space £^(M^, d^x.) (S)£^(M'^, (i^y). Let 
IIo be the projector onto the lowest Landau band, and let n> = 1 — ITo. We will 
use the decomposition 

Hp-(3 = uo{Hp - mo + ^>{Hp - m> 

and 



-^ = no^Ho + n>^n> + no^n> + n>^no. 

X X X X X 



The off-diagonal terms can be bound using 



v^Ho - ^n> ) -^ (^/EUo - ^n> 1 > 

"e / X V ^/e 



for some < e < 1, with the result that 



1 1 / 1\ 1 

— < (1 + £)no^no + 1 + - n>— n>. 

X X V £/ X 



In the same way one shows that, see |[LSY94a 
1 



> (l-3£)no®no- 
3' 



1 



n^®n 



x-y| 
1 



rHn ® Hn 



'^> 



>i rn> 

x-y| 



+ 



e ^ 
l-2£--j Ho^n 

2\ „ „ 1 



W ix> 

^y |x-y 



n^ 



■Hn ® n. 



1-e- 



e 



n-^on 



'^> 



01 



rH^ ® n, 



0- 



(3.5) 
(3.6) 



(3.7) 



(3.8) 



(3.9) 
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Moreover, we will use 
1 



Tro 



n#®n>- 



|x-y| 



n# ® n>r ® r 



< Tr[n#r]supTr 



n>rn>- 



X — a 



(3.10) 

where # stands for either or >. If we restrict ourselves to considering F's with 
Tr[r] < A we therefore have 



£^^[T] > (1 - 3e)^^"[norno] + 36\i-] inf spec \Hp-j3 



2x 



+ Tr[n>r] ( - + inf inf spec {\H0 - -^ 

\2 a \^ £ |x| 

where we have used that 
and that 

^ ^tf3- p- —— I iio ^ Olio I ^/3 



2 1 2A 1 



6 X — a 



Uo{Hs-P- J-] Ho > 5^no f i7« - /3 - i 1 Hn 



2x 



2x 



,(3.11) 

(3.12) 
(3.13) 



for 6 > 1, which can easily be seen by scaling z ^ 6 ^z. 

Now we use the comparison with the hydrogen atom, Proposition (|2.1| ): 



inf spec iHp-/3--—\ > max{l, A"^}E^"(A, /3). (3.14) 

By the same argument as in the proof of this inequality, we can set a = in 
(|3.11|) . Using the diamagnetic inequality we see that |if^ — c|x|^^ > — |c^, so 
finally 



S^m > (l-3E)f™|n„rn„]+3El|l max{l,A}£;™(A./5) 

+ i(,_l(l±^)^.|n,r]. (3,15) 

Now if P is large enough, we can set e^ = 4(1 + XY/P to conclude 
3.2. LEMMA (Comparison with the confined Hartree theory). 

E^^,{X,P)<E^^{X,P){l-const.{l + Xyp-'/'), (3.16) 

where 

E^^„f(X, (3) = 

Tr[nor]<A 



^rf(A,/5) = _mf ^_^,^«[norno] 



(3.17) 
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Note that the constant in ( |3.16| ) can be chosen such that ( p.l6| ) is vahd for all 
values of /? > and A > 0. 

3.3. Remark. Equation ( p.l5| ) can be used to estimate the part of the Hartree 
state r^ not confined to the lowest Landau band. In fact, if e"^ > 4(1 + A)^//3 in 
( |3.15| ), and with F = F^, we get 

Tr[n>FH] < -const.(l + A)^^^^^^^. (3.18) 



Using the simple lower bound on E^^, given in Lemma p.3|, and optimizing over 
e gives 

Tr[n>F"] < const.(l + Xfp-^^^ (3.19) 

3.3 Confinement for the mean field Hamiltonian 

An analogous result as in the previous subsection holds also for the linearized 
theory. The following estimate will be used in Subsection [4.4|: 



3.4. LEMMA (Confinement for the mean field Hamiltonian). Let H be 

given as in l\2. 54\ ), with p > and f p < X. Then for all P > 



inf spec H > inf spec UoHUo + (2 + X)(3-^/^E^^'^{p, 2). (3.20) 



Proof. For < e < 1 we use again ( |3.5| ) and (PT^), and the analogous inequality 
in the other direction for |x|~^ * p, to conclude that 

H > Uo(Hp-P-{l + e)^^+{l-e)^*p]uo 

V m m J 

+n> (hp-p-{i + e'^)^ + (1 - E-^)^ * p) n> 

> (i-£)no//no + £no('ff/3-/3-^jno 

+ |n> [Hp + (3-— — * p n>, (3.21) 



where we have also used ( p. 121 ). With the aid of the inequality ( ^.3] ) one easily 
sees that 

inf spec (h3-%--%-*p\ > inf spec ( Ha - (i±^^)^^ . (3,22) 

V |x| |x| J \ |x| y 

Using again ^if^j — c|x|^^>— |c^we finally get 

H>{1- e)IioHIiQ + eE''^\p, 2) + in> (/5 - (2 + A)^^-^) . (3.23) 

In particular, if we choose e = (2 + A)/9^^/^, we arrive at the desired result, as 
long as (3 is large enough to ensure e < \. But if /3 < (2 + A)^, ( p.20| ) holds 



trivially because of the positivity of p. D 
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3.4 The limit of very strong magnetic fields 

In the investigation of the hmit /? ^ oo some special concepts from earher studies 
are involved: The right scaling, the restriction to the lowest Landau band, the 
effects of the superharmonicity on the distribution of the density in the perpen- 
dicular directions, and the high field limit of the Coulomb interaction. 

We consider pure states A = |x) (xl defined by the wave functions with L" = 
in the lowest Landau band. 



XW 




'/3(x^)' 



''/^Vl^{lJ 



(3.24) 



We denote here the perpendicular and parallel components of the coordinates as 
X = (x-'-,a;"). The scaling factor L = L(P) has been defined in 
solution of the equation 

f3^/^ = L{(3)smh[L{f3)/2]. 

Note that 

L{(3) = In (3 + 0{\n\n (3) 

cxD. In the following we use the scaled coordinates 

z = L{l3)x\ r = a/px"*", r = |r|. 



as /5 



(3.25) 
(3.26) 

(3.27) 

With the density matrix A = |x)(x|) we calculate the contributions to £^^^[A]. 
The normalization of A, we denote it as A^, is equal to HV^Hl- We restrict the set 
of ip to real valued wave functions, since they minimize the kinetic energy, if the 
density \ip\'^ is kept fixed. 

3.5. LEMMA (The strong field limits). Given the state A defined by real 
valued ip G 7^^(R) as in ( 3.24 ), the contributions to the energy in the mag- 
netic Hartree theory, magnetic-kinetic energy Kp^^, attraction Ap^^, and repulsion 
Rl3,ipj fl'^e in the following way related to the kinetic energy, K^ = J{dip/dzydz, 
attraction- energy A^ = ipi^Y > '^^^ repulsion- energy R^ = ^ J il!{z)'^dz of the 
density ipiz)"^ in the HS-theory: 





Km 


= L'K, 


2^2 "^AV- 


-A^ 


^?(^ 


^2 ^Z^''/' 


— R^ 


- L 



1p1 



1/4^3/4 



A, ^ \^/^K^/^ 



(3.28) 
(3.29) 

(3.30) 



Proof. The equation for the kinetic energy is obvious by definition of A. We 
calculate the energy of attraction as 



A 



M 



L' 



-rV2 



Vi3^r{z)^'^{z)dz 



rdr. 



(3.31) 
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where 



V^,r{z) 



For the term in square brackets we use Lemma 2.1 of [pSYOO|| , 

I [...]- ^(0)^1 <A,/r + 8A;/VV^ 
to estimate the difference to the expected hmit in the HS-theory: 



— y4/3,^ - A^ 



For the repulsion Rs ^ we calculate 



L2 



Ri3,ip ~ Rip 



< 



^27r^ ^ 



where 



1 



< 



L 



^ + 8XfKf\r-r'\y' 



r — r 



(3.32) 



(3.33) 



< i (a^ v^ + 8 ■ 2i/^r(5/4)Af irf ) . (3.34) 



(3.35) 

(3.36) 
(3.37) 



The inequality is supplied by Lemma 2.2 of [[BSYOOl] , with an adaptation due 
to the different notation concerning the normalization. After inserting ( p.37| ) in 
( |3.35|) , the integrals can be evaluated as 



1 
Z2 



Rp,ip ~ Ri' 



<^((v^/4)Aj + 4v^r(5/4)-Afir, 



5/4 j^3/4 



(3.38) 
D 



3.6. THEOREM (Magnetic Hartree energy for large /3). In the limit 
/5 ^ oo with A fixed, 

Proof. Combining all the bounds of the Lemma, and using Xj KJ < ^X^ + ^K^ 
to simplify, gives 



1 
L2 



-^ri^]-^""'"""^ 



< 



C 



HP) 



[l + X^){X^ + K^). 



(3.40) 



For the upper bound on E^^{X,l3) we specify ip as ip^^, the minimizing 
wave function for the HS-theory. See equation (3.6) in [ [LSY94a[ ] or the following 
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( ^.30| ). It remains 'ip^^y fo^' all A > 2, so these wave functions are normalized 
^'^ llV^lli = -^ip = min{A,2}. Since the variational principle with fixed norm can 
be replaced by a variational principle with bounded norm, as we have remarked 
in the Subsection |2.1| , the wave function ip^^ is good for the upper bounds for 



all A > 2. K^ is the kinetic energy in the hyperstrong theory. Using the virial 
equation ( |2.42|) , we may replace K^ on the right side of ( |3.40|) by |i?^^(A)|. With 
this choice of ip the equation ( |3.40D gives the upper bound 



-iEnA,/3)<i?«^(A) + ^(A+|E»^(A)|). 



(3.41) 



To derive a lower bound, we use the error bound, when confining the theory 



to the lowest Landau band, which we have estimated in the Subsection 



Lemma 3.2 we have, for /? large enough. 



E^^iX,P) > (1 - const.(l + \yp-'/')-'EZ,i\,P). 



By 



(3.42) 



Also this confined Hartree theory has rotation invariant minimizers, since 
the lowest Landau band is mapped onto itself by rotations around the z-axis. 
Moreover, since the potential is superharmonic, also the minimizer of this confined 
theory has L II = 0. See Corollary ^TTT| of Subsection pl5| . The variational principle 
can therefore be restricted to states A = |x)(x|) where the wave functions x ^i-re 
specified as in equation ( |3.24| ), with general ip G 7i^(R), normalized to ||^||2 = A: 



Konfi^^P) 



inf S^^ 

A, Tr[A]=A, A=|x>(xl ^ 



[A]. 



(3.43) 



The comparison with HS-theory in ( ^.40 ) is now used as a lower bound 



HSr 



HP) 



£r[M>£''v\ 



c 



2^P 



Li(3) 



'1 + \){X + K^). 



(3.44) 



The right side of this inequality can be considered as a functional similar to the 
HS-functional, but with the constant 1 — C(l + A)/L(/3) in front of the kinetic 
energy. With an appropriate scale transformation, this is equivalent to a HS- 
functional multiplied with (1 — (1 + X)C / L{(3))~ , as long as this parameter is 
positive. Taking the infima of both sides, we conclude that 



<nf(A,/5)> 1-C7 



,(1 + A) 



-1 



for (3 large enough. Considering the limit (3 

1 



^"^(A) - C 



(1 + A)A 



L(/3) 
oo at constant A we infer 



"rs^ii^^™ 



(A,,3)>B"'(A). 



Combining this with ( |3.42D proves, in union with ( |3.41| ), the theorem. 



(3.45) 

(3.46) 
D 
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< 1-C— — E"S(A) + C— — , (3.47) 



3.7. Remark. To be precise, we specify the asymptotics: 

(3.48) 
as long as the terms in braces are positive. 

3.8. Remark. The convergence both of the energy of the atom, and of the 
energy per unit charge, when divided by (ln/3)^, is uniform in A on bounded sets 
of A. 

Analyzing the limit A^ ^ oo of many particle quantum mechanics, we will be 
led to the linearized mean field theory, where we need the 

3.9. LEMMA (Generalized strong field limits). Given two states which are 
defined by real valued ip G 7i"'^(M) and real valued ip G H^(M) as in ( ^.24) , with 



densities p/3^^ = P/27ie~'^^^ ^ ^'^Lip{Lx^^y, the following estimate for the Coulomb 
repulsion holds: 



Proof. As for Rp^ in the proof of Lemma |3.5| , with ipi^)'^ ^^ ( p.36|) changed to 

4 The mean field limit 

We now prove the theorems that we have stated in the introduction. We will 
derive appropriate upper and lower bounds to the quantum mechanical energy of 
symmetric states obeying Bose statistics. 

The atom with many interacting particles will be compared to models with 
independent particles in an effective field. In this comparison the upper bound 
is rather easy: One uses symmetric wave functions of product form as trial wave 
functions. Then one adds the "self energies" of the one-particle densities. 

The production of lower bounds is not so easy. For the Regions 1 and 2 we will 
use the Lieb-Oxford bound, [[L081|l , and then we will borrow a part of the kinetic 



energy to estimate the correction in comparison to the mean field model. But 
this method breaks down in the presence of magnetic fields, when they are too 
strong. In Region 3 its use is restricted to the subregion (3 < N'^^^. This situation 



is similar to the case of fermionic electrons, discussed in Sect. 7 of ||LSY94a 



For P > N^^^ we have to borrow some kinetic energy at the level of quantum 
mechanics for A^ particles. Here we develop a new way of producing bounds, 
extending the method of | |BSYOO |. 
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4.1 Upper bounds 

4.1. LEMMA (Upper bounds for Regions 1 and 2). The Hartree energies 
provide upper bounds to the quantum mechanical energies for each N, A and j3: 

E{N,X,(3) E^^{X,(3) 

N - X ' ^ ' ' 

Proof. We apply the variational principle, using A^-fold products of 1-particle 
states r as A^-particle states 

E{N, X,(3) < NTT[{Hp - I3)T] - N J ^^pr + X^^^^^D[pr, pr] (4.2) 

for all r with Tr[r] < 1. Setting F = F^/A we see that the inequality holds, for 
each P and A. D 

4.2. LEMMA (Upper bound for Region 3). In the limit (3 — »• cxd, the en- 
ergies of the hyperstrong theory are asymptotic upper bounds to the quantum 
mechanical energies for each N: 

E{N,X,(3) ^ E^^X) 
lim sup < r . 4.3 

/3^oo N{\n/3y X 



To be precise: 



^W^.«<fl_cJW^^ + ^. (4.4) 



NL{/3)^ - V HP) J X L{l3) 



Proof. This results from combining Lemma 4.1 with Theorem p.6|. D 



4.2 Lower bounds for Regions 1 and 2 

4.3. LEMMA (Lower bounds for Regions 1 and 2). In the limit N —^ cxd, 
the Hartree energies are asymptotic lower bounds to the quantum mechanical en- 
ergies, 

,ta,„f^WM)>^!!!(M), (4,5) 

N->oo N - X ^ ^ 

uniformly in (3 for bounded (3. 

Proof. We use the Lieb-Oxford inequality ||L08! 



^1 Y, Ix* - x.r'l^) > D[p^,p^] - C I pf (4.6) 
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(\1/ is normalized as ||^I/||2 = 1, the constant can be chosen to be C = 1.68), which, 
together with Holder's inequality J p'^^^ < {j P^Y^^{j pY^^i implies that 

(v|>|i7^,,,^vi>) > ^Ef''[T^\/N] - CXN-'/' (jpl\ , (4.7) 

where V^, is the one-particle reduced density matrix of \l/, and p^ is its density. 
Now if (\E'|ifiv A,/3\E') < 0, which we can of course assume, then 

{M Y.^Hp, - /?)vl>) < -{^ Y.^Hp,, -a- 2|x,r')^) < ~NE'^^\(i, 2), (4.8) 



with E^y'^ defined in (^^. Together with (|2:20| ) this implies that 

1/3 1 /o\4/3 



(/p|) ^\{^ N{(3-E^y^{(3,2)), (4.9) 



so finally 



^E{N,X,(3) > ^E'^'^iX, p) - CXN-^/^p - E''^^{(3,2)Y^\ (4.10) 

A* A 



D 



4.4. Remark. Note that the convergence of the energies in Region 2, including 
Region 1, is uniform in (3 for bounded jS. So if /5 — * as A^ ^ oo we get the 
usual Hartree energy without magnetic field. It is even possible to let /3 ^ cxd 
with A^ as long as /5 < const. N^^^. In fact, in ( [4.10|) we have an error term of 
order A^~^/^/?^/^ (note that E^^'^ ~ (In/?)^ for large /3), and this is of lower order 
than E^^ as long as /3 < const. iV^/^. 

4.3 Restriction to independent particles 

An essential ingredient is the positive definiteness of the repulsive pair interaction. 
The method exploiting positive definiteness of a function W uses the inequality 

^iy(x,-x,.) > J2 W^(x, - y)a(y)dV (4-11) 



Actually we need a function W which is positive definite, finite at the origin, 
and a lower bound to the Coulombic repulsive potential. Also the cutoff near the 
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origin should vanish in the hmits A^ -^ oo, f3 — > oo, but still with W{0)/N{\n(3)'^ 
going to 0. If we choose for W the spherical symmetric cutoff potential 

I _ g-M|x| 

Kutoff(x) = rn , (4.12) 

it turns out that this does not suffice for our purposes. The cutoff length /i~^ 
should be of the order of the typical lengths of the atom, and this would require 
a coupling of the limits N ^ oo and /? — i> oo, or more precisely, (3 is not allowed 
to increase arbitrarily fast with A^. 

To get a useful lower bound for the entire Region 3, we have to push down 
ly (0) even further, without changing the effective interaction to much. To achieve 
this, we split Kutofr (we restrict the parameter to fi > 1) into its effective part V^ 
and the long range tail Viong, 

e-|x| _ p-M|x| 
^m(x) = rn , /">1, (4.13) 

Mong(x) = ^—^, (4.14) 

and borrow some kinetic energy. 

The construction of W proceeds in several steps. We begin with the ob- 
servation that V^ is integrable in x" for any x^. It has moreover the property of 
being decreasing in |x|. So the integrals along lines with fixed x-*- can be bounded 
by the integral along the line, where x-*- = 0: 

/oo 
V;.(x)t/xll <21n/i. (4.15) 

'OO 



We now use Lemma 6.3 of [PSY00|| , the operator inequality, which holds 



for a > 0, F > 0, and for each b > 0, b independent of a and F, 

apl + FSix) > F ■ WF/a,b{x), (4.16) 

where 

wn,,{x) = D^e-'^\^\ (4.17) 

and p"^ = —{d/dxY. These functions have the properties of being positive definite 
and producing delta-function sequences in the limit bD — > oo coupled with b — > 
oo. (The variable x will in the application be replaced by x".) We extend this 
lemma to the 

4.5. PROPOSITION (Operator inequalities). Let V(x) define a positive 
Borel measure V{x)dx, with j_ V{x)dx < F, F > 0, a > 0. For each b > 
the operator inequality 

apl + V{x) >{V* WF/a,b)ix) (4.18) 
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holds, in the sense of an inequality for quadratic forms, with the Sobolev space 
7i^(M) as the form domain. 

Proof. For ip G 7i^{M.) consider the wave functions ipy{x) = ipi^ + y). By ( |4.16| ) 
there holds for each y the inequahty 

{^y\{ap^ + F6)\^y) > F{^y\wF/a,b\^y). (4.19) 

Both sides of this inequahty are Borel measurable functions of y. Observe 

(ipylp'^li'y) = ii'lp'^li'), {i^y\w{x)\ipy) = (V' | w (x — ?/) |'?/') , sud lutcgrate both sides 
of ( |4.19| ) with V{y)dy. Divide by F, and use j_ °° V{x)dx/F < 1 in front of the 
kinetic energy term. D 

We may consider a{plY + V"^(xj — Xj) as a 1-particle Hamiltonian acting 
in C'^(R,dxl), which is parametrized by the perpendicular coordinates and by 
Xj. However, the operator inequality ( |4.18| ) is also valid if the operators act in 
the extended space £^(M'^^, (i^xi...(i^X7v), so we can apply this proposition, with 
F = 21n/i, to 

a(pf)2 + \/^(x,-x,). (4.20) 

We add the missing long range tail and define 

W{^) = [ V,{^ - y)S\y^)w2in^./a,b{y^^)d'y + ^V"' (^'^l) 

J P^ 



where w,„{^ ) is defined in ( [4.171) . The essential effect of the convolution of V^ 
with the distribution 

w(x) = (5(x^)w;(x!l) (4.22) 

is the lowering of the value of the interactions at the points of coincidences of 
particles, 

W{0) < (jv^dx\\\ ^21n;./a,6(0) + 1 < ^(l^/x)' + 1- (4.23) 

Thus we have shown, that W^(x) is a lower bound to the Coulomb potential, 
when some borrowed part of the kinetic energy is added, 

W = V^*w + V,,^^ < a{p\\f + V^ + Viong < a{p\\f + l/|x|. (4.24) 

Now the function W = V^*w + Vio^g is also positive definite, since the distribution 
w and all the involved functions are positive definite. So we can refer to the 
technique ( |4.11D of reduction to one-particle models. 
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4.6. PROPOSITION (Bounds by independent particles). Considered as 
quadratic forms, the N-particle Hamiltonians Hn,\,i3 are bounded from below by 
sums of one-particle Hamiltonians and constants: 

Hnxp>Y.^^-~ f f cr{^)W{^-y)a{y)d'^d'y - X-{\nf,f - ^, (4.25) 



where 



h, = Hp, - /5 - '^{plf - i^ + A I w{^, - y)a(y)ciV, (4.26) 



and where W^(x) is defined in ^.21\) . The parameters a, h have to be positive, 
fi> 1, 0"(x) should be an element of C^(M.^). 



Proof. From Hf^^i we borrow a part of (p|)^ = —{d/dx-)"^ to use ( [4 .241 ), 

a(pl')' + , ^^1 > W^(x. - X,), (4.27) 

sum over all pairs i y^ j , and multiply with X/2N. This gives 



Inserting (|4.28|) and ( [4.11|) into (|1.3| ) completes the proof. D 



To ascertain any use to the inequality, aA < 2 is requested. In the next sub- 
section, the density a will be chosen to be an approximately minimizing Hartree 
density, and the parameters will be adapted to the limit N, j3 ^ oo. 

4.4 The lower bound for Region 3 

We choose a = Np, with 

p(x) = Aexp (-^) \LpTiLx^ (4.29) 



The minimizing density p^^ for the hyperstrong theory is, see |[LSY94a 



tanhc(A) = (2 - A)/2, 
p^^{z) = 2{2+\z\)-^ for A>2. (4.31) 
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Note that ||p||i = minjl, 2/A}, ||p||oo = (/9/27r)L(/?)pP(0)/A < C(3L{(3). 

We want to study the lower bound to N^^E{N, A, P), due to Proposition |]6, 



inf spec {h} - ^ jj p{-^)W{-K - y)p{y)S^d'y - ^ (^(In/i)' + ^) , (4.32) 

which we expect to be asymptotically proportional to (ln/5)^. For simplicity, we 
will restrict our considerations io (3 > Cp for some Cp > 1 in the following. We 
choose the parameters as 

/, = /3i/2+'5, a = A^-^ 6 = AT'', (4.33) 

with (5, £, 7] all greater than zero, e + ry < 1. This guarantees first of all the relative 
vanishing of the constant which stems from iy(0): 

^(^(ln/i)^ + 0/(ln/3)^-O, (4.34) 

as f3,N —* oo. It remains to study the Hamiltonian h and the self energy of p 
due to the interaction W = V^* w + V^ong- 
The inequality 

jj p{^)iV,*w)i^-y)p{y)d'^£y < j j p{^)V,{^- y)p{y)d?^Sy (4.35) 

can be observed in Fourier space, where V^ and w_ are positive, and w_<l. We add 
l^ong, and use the pointwise inequality in x-space (V^+Viong)(x) < V^{'k) = l/|x|. 
So the self-energy of p due to W is smaller than the self-energy D[p, p] due to the 
interaction by the Coulomb potential V'-'' . 

The potential (V^ * w + V[ong) * p in the Hamiltonian h is now changed to 
V" * p. First we give a bound to the difference of V^ * w * p to V^ * p in the form 
of an operator inequality, which follows from 

\{ij\Vf,*{w*p-p)\ij)\ < ||(|^|^*v;)-(w*p-p)||i 

< \\W* p- p\\i II 1^1^*1^^1100 

= ll^*p-p||i sup(^|V^,y|V^), (4.36) 

y 

where V^,,y = V;,(x - y). Since V^ < V^, and, with E^y<^(/?) = E^'y^{f3, ( = 1), 

Hp-l3-V^>E''y\[3), (4.37) 

the last term in ( ^4.36| ) is bounded by 

{^\V^\^) <{mp-(3- E^^^^if^m, (4.38) 
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multiplied with \\w*p — p\\i. 



We observe, see Remark 2.2 



^^ ^ " A A^o A 4 



A-+0 



The precise bound in the Remark ^]^ gives thus (see also ||AHS81|| ) 



(4.39) 



(4.40) 



for P > Cp. An analogous bound holds also for E^^'^{P,2). Note that L{P) can 
be replaced by (In/?) (and vice versa), since for /5 > C/j > 1 there are positive 
constants Ci and C2 such that Ci < ln(/3)/L(/9) < C2. 

To complete the estimate ( [4.36| ), we need a bound to ||w*p — p||i. There, the 
integrals in perpendicular coordinates can be done explicitly. Note that for all A 
we have \dp^^/dz\ < p^^. Using this and the monotonicity of p^^ in \z\ we can 
estimate 



\pxiz-y)-pxiz)\ < \y\ sup 

xe[z-y,z] 



d^ 

dz 



[x] 



< \y\ [pT{z -y)+ pT{z) + Q{\y\ - k|)pf (0)) . 

(4.41) 



Since for \y\ > 1 ( [4.41|) is obviously true even without the last term, we see that 
(^4.41|) holds with Q{\y\ — \z\) replaced by 9(1 — |2;|). This implies 



1 /• II 
||p*ii; — p||i < — dx"driLw{ri) 

X (|pf (La:" - Lry) - pf (La:")] + {{J w)-' - l) pf (Lx")) 

< A- V pf (1 -Jw) + 2A-1 (/ pf + pf (0)) L J \v\wiv)dv. 

(4.42) 

Now /p^s < x^ pHS(o) < lA, {1- Jw) < (26)-i and J \r]\w < a/{2b\np), so we 
get 



\p * "w — p\\i 



< 



1 3 aL{p) 
26 ^2 61np 



< CiN-'^ + N-'^-' L{f3)/ In p) < C N~" . (4.43) 



Finally, we add to V^ the short-range term 



Kihort (x) 



-^l\yL\ 



(4.44) 
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and we find the pointwise bound to (V*" — V^) * p as 

llVshortllillplloo < Cfi-^PL{P) = Cp-^'L{P). (4.45) 

We collect the bounds we have obtained so far: 
E(iVA/3) ^.^f gp^^ {h}-\D[p,p]-C\i\nPf {N'^'^-' + N~^ + p-^'L{P)-') , 

(4.46) 
with the one-particle Haniiltonian 

hh = {Hp - /3) (1 - CAA^^" - 2\N~') -V^ + \V^ * p, (4.47) 

where we have used (p")^ < Hp — (3. We now use ( |4.8|) to estimate the error 
terms in the kinetic energy. Then we apply the statements of the subsection on 



the confinement in the mean field theory, Lemma |3.4| , to estimate 

^^IIjMI > infspec {no(i7;3-/5-y^ + Al-^*p)no}-AD[p,p] 

(4.48) 
Corollary |2.11| states that the minimizer of the operator in question has L" = 0. 



The search for the infimum of the spectrum can therefore be restricted to the 
use of the same type of wave-functions as in the investigations on the limit of 



very strong magnetic fields in Subsection |3.4| , Lemma |3]^ and |3.9| , and these 
investigations can be applied, too. We choose e = t] = 1/3 and get the lower 
bound 

E{N,\,P) 

> (inf spec {pl-6{z)+pT{z)}-^J{pT{z)rdz^ ~ ^iw) ^ 

- C (AiV-i/3 ^ ^pyi ^^^-25 ^i^x) + {2 + A)/3-i/2) . (4.49) 

The operator in curly brackets is the Hamiltonian for the linearized HS theory, 
with ground state energy (E'^^H- ^ /(Pa^)^)/A (cf. eq. ( p.38| )). So our final result 
can be stated as 

4.7. LEMMA (Lower bound for Region 3). If L{P) > C(l + A), then 

NL-m - A y ^m) ^v"" ^ m)- ^'-'^^ 

Hence the convergence of the lower bound, in the limit N —^ oo, (3 -^ oo, is 
proven. 

4.8. Remark. The limits N ^ oo and /3 — ^ cx) can be considered independently 
of each other, in any succession or combination. 
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4.5 Critical particle number 

Recall that E_[N, Z, B) denotes the ground state energy of the unsealed Hamil- 



tonian ([LI 



4.9. LEMMA (Limit of the critical charge). Define 

N^{Z, B) = max{N\E{N, Z, B) < E{N -1,Z,B)}. 
We have 

Z^oo Z ~ 

Proof. This follows from the convergence of the energies, by analogous arguments 
as in [[BL8l or [|LSY94a|| . 



(4.51) 



(4.52) 



D 



4.10. Remark. In [pOCI|] , Theorem 3, the following upper bound on Nc is proven: 



A'c<2Z + l + |mmi (l + ^|,C|l 



In (4 

\Z^ 



(4.53) 



for some constant C independent of B and Z. Inserting ([4.53|) into ( [4.52 ) this 



proves Lemma 2.£ 



4.6 Convergence of the density matrices 

4.11. LEMMA (Ground states are Hartree minimizers). For fixed X and 
(5, let '^ N he an e-approximate ground state of Hn,x,i3 defined in dl.dj ), i.e. 



(^;v, ^^,A,/3^7v) < EiN, A, P) + eiN) 



(4.54) 



for all N, with < e = o{N). Let r^v be the corresponding one-particle reduced 
density matrix. Then XTj^j/N is a minimizing sequence for S^^, i.e. 



lim Sf'^lXTN/N] = E^^(A,/3). 

Proof. Let pn be the density of T^. We have 

A 



(4.55) 



A2 



^MH < S^n^XTr^/N] = -Tr[(^, - /? - |x|-^)^^.] + ir2D[PN,PN 



N 



m' 



X x^ x^ 

-A^>n\Hnxp'^n) - T7^(^^l Y. 1^* - ^il^'^A^) + JrlD[pN,PN] 



N 



< ^{E{N,X,(3)+e) + ^^C I p% 



i<j 

4/3 



(4.56) 



where we used again the Lieb- Oxford inequality for the last step. By an analogous 
argument as in Lemma ^^ we see that the right side of ([4.56[ ) converges to E^^ 

as N —^ oo. n 
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Since XT^/N is a minimizing sequence for S^^, and F^ is unique, we know 
that XTn/N -^ r^ in weak operator sense. If A < Ac, there is even norm 
convergence. To show this we need the following general lemma: 

4.12. LEMMA (Weak convergence implies norm convergence). Let 

aN he a sequence of positive trace class operators on a Hilbert space. Suppose 
that aj^ ^ a in weak operator sense, for some positive trace class operator a. If 
linijv^oo Tr [fljv] = Tr[a], then a^ ^ a in trace norm, i.e. 

II^Af — a||i = Tr[|aAr — a|] — > as A^ -^ oo. (4-57) 

Proof See Em, Thm. 2.20. D 



4.13. THEOREM (Convergence of the density matrices). Let Tn be as 



in Lemma [^ . i 4 ^'^'^ l^i A < Ac- Then for each fixed j3 

lim ||Ar^/A^-r^||i = 0. (4.58) 

Proof. Since F^ is unique, XT j^/N converges weakly to F^ by Lemma [4.11| and 



the proof of Theorem |2.4| . Apply Lemma [4.12| to a^ = XT^/N and a = F^, and 
note that Tr[AFjv/A] = A and Tr[F"] = min{A, AJ. D 

4.7 Bose condensation 



Bose condensation, as defined in [ PU5(j|| , means that the largest eigenvalue of the 



one-particle reduced density matrix of the ground state for the A-particle Hamil- 
tonian is of order A as A — > cx), or equivalently, if liminfTv-^oo l|rAr||/Tr[F7v] > 0. 
This is shown to be the case for our model. 

4.14. THEOREM (Bose condensation in the mean field limit). Let T^ 

be as in Lemma \4 . 1 4 Then for each fixed (3 and X 

lim JIM = '^^^. (4.59) 

^^ooTr[Fjv] A ^ ' 

In particular, if X < Ac, and ififN denotes the normalized eigenvector correspond- 
ing to the largest eigenvalue of Fjv with appropriate phase-factor, we have 

IIP II 
lim .^ , = 1, lim ||(^7v-V^h||2 = 0, (4.60) 

N^oo lr[i TvJ N^oo 

where (fu = X~^^'^^/p^ is the normalized ground state of H^. 
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Proof. We use the notations of the proof of Theorem |4.13| . The first assertion is 
easily proved, using that 

|||a7v|| — ||a||| < llctAf — a\\ < \\aN ~ fl||i (4.61) 

and ||a|| = min{A, Ac}. To prove the second we denote P = |(^h)(v^h| and P/v = 

\'^n){'^n\, and compute 

Tt[aNP] = \\aN\\{^N\PVN)+TT[a]i^Pa]i\l-PN)] 

< ||aAr||(v3iv|P(/?Ar) +Tr[aAr] - llaivll, (4.62) 

where we have used P < 1 in the last step. Therefore {(Pn\P^Pn) ^ I as N ^ oo, 
which gives immediately the desired result. D 
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